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00 ' We study constraints on additional Z' bosons predicted in the supersym- 

^ ■ metric (SUSY) models by using the updated results of electroweak ex- 

p • periments - Z-pole experiments, my/ measurements and low-energy neu- 

p tral current (LENC) experiments. We find that the effects of Z-Z' mixing 

^ ■ are parametrized by (i) a tree-level contribution to the T-parameter, (ii) 
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the effective Z-Z' mass mixing angle ^. In addition, the effect of the di- 
rect exchange of the heavier mass eigenstate Z2 in the LENC processes is 
parametrized by (iii) a contact term g'Elc^'m?z2- S^^^ theoretical pre- 
dictions for the observables in the electroweak experiments together with 
the standard model radiative corrections. Constraints on Tnew and from 
the Z-pole and mw experiments and those on gE/c^rn"^^ from the LENC 
experiments are separately shown. Impacts of the kinetic mixing between 
the U(l)y and U(l)' gauge bosons on the x^-analysis are studied. We show 
the 95% CL lower mass limit of Z2 as a function of the effective Z-Z' mixing 
parameter (, a combination of the mass and kinetic mixings. Theoretical 
prediction on ( and Qe is found for the Xii^iV aiid z/ models by assuming 
the minimal particle content of the SUSY Eq models. In a certain region of 
the parameter space, the Z2 boson mass in the detectable range of LHC is 
still allowed. 
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1 Introduction 



Although the minimal Standard Model (SM) agrees well with current electroweak 
experiments it is important to examine consequences of new physics models 
beyond the SM at current or future collider experiments. One of the simplest 
extensions of the SM is to introduce an additional U(l) gauge symmetry, U(l)', 
whose breaking scale is close to the electroweak scale. The U(l)' symmetry is 
predicted in a certain class of grand unified theories (GUTs) with gauge group 
whose rank is higher than that of the SM. In general, the additional U(l)' gauge 
boson Z' can mix with the hypercharge U(l)y gauge boson through the kinetic 
term at above the electroweak scale, and also it can mix with the SM Z boson 
after the electroweak symmetry is spontaneously broken. Through those mixings, 
the Z' boson can affect the electroweak observables at the Z-pole and the W boson 
mass mvi/- Both the Z-Z' mixing and the direct Z' contribution can affect neutral 
current experiments off the Z-pole. The presence of an additional Z' boson can 
be explored directly at collider experiments. 

The supersymmetric (SUSY) models are the promising candidates which 
predict an additional Z' boson at the weak scale (for a review, see 0). The gauge 
group i?6 can arise from the perturbative heterotic string theory as a consequence 
of its compactification. In the models, the SM matter fields in each generation 
are embedded into its fundamental representation 27 that also contains several 
exotic matter fields - two SM singlets, a pair of weak doublets and color triplets. 
Because E^ is a rank-six group, it can have two extra U(l) factors besides the 
SM gauge group. A superposition of the two extra U(l) groups may survive as 
the U(l)' gauge symmetry at the GUT scale. The U(l)' symmetry may break 
spontaneously at the weak scale through the radiative corrections to the mass 
term of the SM singlet scalar field 0. 

In this paper, we study constraints on the Z' bosons predicted in the SUSY E^ 
models. Although there are several previous works 0, |^, ||, |^, |[ , we would like to 
update their studies by using the recent results of electroweak experiments, and 
by allowing for an arbitrary kinetic mixing |10|, [Tl|] between the Z' boson and 
the hypercharge B boson. In our study, we use the results of Z-pole experiments 
at LEPl and SLC, and the m-w measurements at Tevatron and LEP2 which were 
reported at the summer conferences in 1997 We also study the constraints 
from low-energy neutral current (LENC) experiments: lepton-quark, lepton-lepton 
scattering experiments and atomic parity violation measurements. 
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We find that the lower mass hmit of the heavier mass eigenstate is obtained 
as a function of the effective Z-Z' mixing term which is a combination of the 
mass and kinetic mixings. In principle, C is calculable, together with the gauge 
coupling once the particle spectrum of the model is specified. We show 
the theoretical prediction for ^ and in the SUSY models by assuming the 
minimal particle content which satisfies the anomaly free condition and the gauge 
coupling unification. For those models, the electroweak data give stringent lower 
mass bound on the Z^ boson. 

This paper is organized as follows. In the next section, we briefly review the 
additional Z' boson in the SUSY E^ models and the generic feature of Z-Z' mixing 
in order to fix our notation. We show that the effects of Z-Z' mixing and direct Z' 
boson contribution are parametrized by the following three terms: (i) a tree-level 
contribution to the T parameter Tnew, (ii) the effective Z-Z' mass mixing angle 



^ and (iii) a contact term gE/Cx^'z2 "which appears in the low-energy processes. 
In Sec. 3, we collect the latest results of electroweak experiments. There, the 
theoretical predictions for the electroweak observables are shown together with 
the SM radiative corrections. In Sec. 4, we show constraints on the Z' bosons 
from the electroweak data. The presence of non-zero kinetic mixing between the 
U(l)y and U(l)' gauge bosons modifies the couplings between the Z' boson and the 
SM fermions. We discuss impacts of the kinetic mixing term on the x^-^nalysis. 
The 95% CL lower mass limit of the heavier mass eigenstate Z2 is given as a 
function of the effective Z-Z' mixing parameter (. The ^-independent constraints 
from the low-energy experiments and those from the direct search experiments at 
Tevatron are also discussed. In Sec. 5, we find the theoretical prediction for ( 
in some SUSY Eq models (x, ip, i^) by assuming the minimal particle content. 
Stringent Z2 boson mass bounds are found for most models. Sec. 6 summarizes 
our findings. 

2 Z-Z' mixing in supersymmetric Eq model 
2.1 Z' boson in supersymmetric Eq model 

Since the rank of E^ is six, it has two U(l) factors besides the SM gauge group 
which arise from the following decompositions: 

E, dSO(10)xU(1)^ 

D SU(5) X U(l)^ X U(l)^. 
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An additional Z' boson in the electroweak scale can be parametrized as a linear 
combination of the U(l)^ gauge boson Z^ and the U(l);^ gauge boson Z^ as [|I3] 

Z' = Z^ cos Pe + sin (3e- 

In this paper, we study the following four Z' models in some detail: 



(2.2) 



I3e 





Tl/2 


tan-i(-^5/3) 


tan~i(\/l5) 


model 


X 




7] 


V 



(2.3) 



In the SUSY-£'e models, each generation of the SM quarks and leptons is embed- 
ded into a 27 representation. In Table 1, we show all the matter fields contained 
in a 27 and their classification in SO(IO) and SU(5). The U(l)' charge assignment 
on the matter fields for each model is also given in the same table. The nor- 
malization of the U(l)' charge follows that of the hypercharge. Besides the SM 
quarks and leptons, there are two SM singlets v'^ and S", a pair of weak doublets 
Hu and if^, a pair of color triplets D and D in each generation. The 77-model 
arises when breaks into a rank-5 group directly in a specific compactification 
of the heterotic string theory |I4[. In the i/-model, the right-handed neutrinos v'^ 
are gauge singlet |T^ and can have large Majorana masses to realize the see-saw 
mechanism ||16|| . 

The U(l)' symmetry breaking occurs if the scalar component of the SM singlet 
field develops the vacuum expectation value (VEV). It can be achieved at near 
the weak scale via radiative corrections to the mass term of the SM singlet scalar 
field. For example, the terms SDD and SHuH^ appear in the SU(3)(^ x SU(2)^ x 
U(l)y X U(l)' invariant superpotential. If the Yukawa couplings of the SDD term 
and/or SHuH^ term are 0(1), the squared mass of the scalar component of 5* can 
become negative at the weak scale through the renormalization group equations 
(RGEs) with an appropriate boundary condition at the GUT scale. Recent studies 
of the radiative U(l)' symmetry breaking can be found, e.g., in ref. 

Several problems may arise in the Eq models from view of low-energy phe- 
nomenology [0. For example, the scalar components of extra colored triplets D 
in 27 could mediate an instant proton decay. It should be forbidden by imposing a 
certain discrete symmetry on the general SU(3)^ x SU(2)^ x U(l)y x U(l)' invari- 
ant superpotential. Except for the //-model [0, the large Majorana mass of v'^ is 
forbidden by the U(l)' gauge symmetry, and the fine-tuning is needed to make the 
Dirac neutrino mass consistent with the observation. Further discussions can be 
found in ref. p|. In the following, we assume that these requirements are satisfied 
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Table 1: The hypercharge Y and the U(l)' charge of all the matter fields in a 
27 for the x, ijj, rj and u models. The classification of the fields in the SO(IO) and 
the SU(5) groups is also shown. The value of U(l)' charge follows the hypercharge 
normali zation. 

SO (10) SU(5) field Y 2V6Q^ ^J72/5Q^ Qr, Qu 




by an unknown mechanism. Moreover we assume that all the super-partners of the 
SM particles and the exotic matters do not affect the radiative corrections to the 
electroweak observables significantly, i.e., they are assumed to be heavy enough to 
decouple from the weak boson mass scale. 

2.2 Phenomenological consequences of Z-Z' mixing 

If the SM Higgs field carries a non-trivial U(l)' charge, its VEV induces the Z- 
Z' mass mixing. On the other hand, the kinetic mixing between the hypercharge 
gauge boson B and the U(l)' gauge boson Z' can occur through the quantum effects 
below the GUT scale. After the electroweak symmetry is broken, the effective 
Lagrangian for the neutral gauge bosons in the SU(2)^ x U(l)y x U(l)' theory is 
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given by |T^ 

+ml^,Z^Z'^ + \mlZ^Z, + ^m|,Z'^Z;, (2.4) 

where F^"{F = Z, Z' , A^ , B) represents the gauge field strength. The Z-Z' mass 
mixing and the kinetic mixing are characterized by and sinx, respectively. 

In this basis, the interaction Lagrangian for the neutral current process is given as 

Cmc = -T.{^ QfJ~al'faK + 9zLr (4 - Qf^ sin^ Ow) LZ^ 

+9EQi"LrfaZ',}, (2.5) 

where gz = g/ cos^vf = dv/ sinOw- The U(l)' gauge coupling constant is denoted 
by gE in the hypercharge normalization. The symbol fa denotes the quarks or 
leptons with the chirality a {a = L or R). The third component of the weak 
isospin, the electric charge and the U(l)' charge of fa are given by J|^, Qf^ and 
Q-^£ , respectively. The U(l)' charge of the quarks and leptons listed in Table 1 
should be read as 



— Qe ~ Qe ^ Qe ~ Qe ~ Qe > \ 
Q^E='Q^i if = e,u,d), /■ 

The mass eigenstates {Zi, Z2, A) is obtained by the following transformation; 



(2.6) 



Z \ I cos ^ + sin ^ sin 6w tan x — sin ^ + cos ^ sin 6w tan x ^ \ ( Zi 

Z' \ = \ sin^/ cosx cos^/ cosx 

A'^ I \ — sin ^ cos 6w tan x — cos ^ cos 6w tan x 1 / \ ^ 



(2.7) 



Here the mixing angle ^ is given by 



tan 2^ = 2 2 2"~\ 2 i o 2 ' (2-8) 

with the short-hand notation, c-^ = cosx, = sinx and sw = sin 9w- The 
physical masses mz^ and mz^ ijnzi < ""^^2) are given as follows; 



= m 



= m|(c5 + sgsvyt^)^ + m|, f — j + 2m|2/ — (cg + sgsvi/tx)) (2.9a) 
lictiSwtx- s^)^ + ml,(^] +2m\z' — {c^swtx-s^), (2.9b) 
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where = cos^, = sin^ and = tan^. The hghter mass eigenstate Zi should 
be identified with the observed Z boson at LEPl or SLC The excellent agreement 
between the current experimental results and the SM predictions at the quantum 
level implies that the mixing angle ^ have to be small. In the limit of small ^, the 
interaction Lagrangians for the processes Z12 fafa are expressed as 



Z2 



JL 



w 



he 



9zC 



w 



9e 



faZ: 



where the effective mixing angle ^ in eq. ( |2.10a|) is given as 

9e 



9z cos X 



(2.10a) 
(2.10b) 

(2.11) 



In eq. ( p.lO|) , the effective U(l)' charge is introduced as a combination of 
and the hypercharge Yf^: 



9e 



(2.12a) 
(2.12b) 



where the hypercharge Yf^ should be read from Table 1 in the same manner with 
Q;^ (sGG) GQ- ( |2.6| )). As a notable example, one can see from Table 1 that the 
effective charge of the leptons (L and e'^) disappears in the //-model if 6 is 
taken to be 1/3 |Tg. 

Now, due to the Z-Z' mixing, the observed Z boson mass m^^ at LEPl or SLC 
is shifted from the SM Z boson mass m^: 

,2 



< 0. 



(2.13) 



The presence of the mass shift affects the T-parameter [0 at tree level. Following 
the notation of ref. |T^, the T-parameter is expressed in terms of the effective form 
factors (7|(0), (/^(O) and the fine structure constant a: 

9wi^) 



aT 



m 



w 



9m 



where Tgj^ and the new physics contribution Tnew are given by: 



1 - 



9w{^) ml 



m 



w 



Am^ 



> 0. 



(2.14a) 
(2.14b) 

(2.15a) 
(2.15b) 
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It is worth noting that the sign of Tncw is always positive. The effects of the Z-Z' 
mixing in the Z-pole experiments have hence been parametrized by the effective 
mixing angle ^ and the positive parameter Tnew 

We note here that we retain the kinetic mixing term 5 as a part of the effec- 
tive Zi coupling Q^^ in eq. ( |2.12aD . As shown in refs. [|Ty, |n|, |TB], the kinetic 
mixing term b can be absorbed into a further redefinition of S and T. Such re- 



parametrization may be useful if the term Y^^b in eq. (|2.12a|) is much larger than 
the Z' charge Q^^ . In the models studied in this paper, we find no merit in 
absorbing the Yfb term because, the remaining Q^^ term is always significant. We 
therefore adopt Q^^ as the effective Z\ couplings and Tnew accounts only for the 
mass shift ( |2.13 ). All physical consequences such as the bounds on ^ and are 
of course independent of our choice of the parametrization. 

The two parameters Tnew and ^ are complicated functions of the parameters 
of the effective Lagrangian ( p.4|) . In the small mixing limit, we find the following 
useful expressions 



i = - 



aTn, 



(9e mz 
\gz rriz' 
Qe mz 
gz mz' 



c 



nir. 



vrir 



l + Oi 



2 n 



where we introduced an effective mixing parameter C, 



c 



9z m^zz' 



-5. 



:2.16a) 



(2.16b) 



(2.17) 



Qe 

The Z-Z' mixing effect disappears at C = 0. Stringent limits on irtz' and hence 
on mz2 can be obtained through the mixing effect if C, is 0(1). We will show in 
Sec. 5 that Q is calculable once the particle spectrum of the model is specified. The 
parameter C, plays an essential role in the analysis of Z' models. 

In the low-energy neutral current processes, effects of the exchange of the heav- 
ier mass eigenstate Z2 can be detected. In the small limit, they constrain the 
contact term g^/c^y^m^z^. 

3 Electroweak observables in the Z' model 



In this section, we give the theoretical predictions for the electroweak observables 
which are used in our analysis. The experimental data of the Z-pole experiments 
and the W boson mass measurement are summarized in Table 2. Those for the 
low-energy experiments [H] are listed in Table 0. 
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^,,11 (data)— best fit 

^ (error) 
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-1.1 
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-1.1 
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1.4 


0.6 
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0.1734 ± 0.0048 
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0.0984 ± 0.0024 
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-2.2 


Ai,c 
^FB 
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-0.3 
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0.0002 


0.0002 


-0.0001 


0.0002 
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Table 2: Summary of electroweak measurements for the Z-pole experiments and the 
mw measurement Q. The best fits to all the data in this Table are found by allowing 
the five parameters m^, as{mzi), Q;(mU, Tnow and ^ to vary freely under the constraints 
rut = 175.6±5.5 GeV H], a/(mzi) = 0.118±0.003 ||], l/a(m|J = 128.75±0.09 |2l|, 
Tncw > and rrijj = 100 GeV. The results for the X) "0) V and u models are obtained by 



setting 5 
5 = 1/3. 



0. The symbol rj* denotes the leptophobic ?7-model where 6 is taken to be 
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^^^11 (data)— best fit 
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175.6 ±5.5 


171.6 
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171.5 
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0.118 ±0.003 


0.1185 


0.1181 


0.1181 


0.1185 


0.1181 


0.1189 
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T 

new 
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128.75 


128.75 
0.0 


128.75 
0.0 


128.73 
0.0 
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0.0 
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0.0 
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0.0002 


-0.0003 


0.0001 
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0.279 


1.771 


-0.646 
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Table 3: Summary of measurements for the low-energy neutral current experi- 
ments [0. The best fits are found by using all the electroweak data of Table 2 
and those in this Table. The results for the Xi'4'iV and v models are obtained by 
setting 5 = 0. The symbol vj* denotes the leptophobic //-model where 5 is taken to 
he 5 = 1/3. 
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3.1 Observables in Z-pole experiments 

The decay amplitude for the process fafa is written as 

^ UZ) = M[ ez, ■ Jf^, (3.1) 

where e^^ is the polarization vector of the Z-^ boson and J^^ = fal^fa is the 
fermion current without the coupling factors. The pseudo-observables of the Z- 
pole experiments are expressed in terms of the real scalar amplitudes with the 
following normalization ||l[] 

,f , Mi Mi (32) 

^AV2G,ml 0.74070- ^ ' ^ 

Following our parametrization of the Z-Z' mixing in eq. ( p.lOa| ), the effective 
coupling in the Z' models can be expressed as 

gi = (gDsM + Q^il (3.3) 

The SM predictions |1^, Q for the effective couplings {qDsm can be parametrized 



as 

{gl)sM= 0.50214 + 0.453 A^l, (3.4a) 

{glhu = -0.26941 - 0.244 Agl + 1.001 As^ (3.4b) 

(^r)sm= 0.23201 + 0.208 A^| + 1.001 As^, (3.4c) 

{gl)sM= 0.34694 + 0.314 A^|- 0.668 As^ (3.4d) 

{g'^)sM = -0.15466 - 0.139 Agl - 0.668 As^ (3.4e) 

(^^)sM = -0.42451 - 0.383 Agl + 0.334 Af, (3.4f) 

{g'k)sM = 0.07732 + 0.069 Agl + 0.334 As^ (3.4g) 

(^l)sm = -0.42109 - 0.383 Agl + 0.334 Af + 0.00043xt, (3.4h) 

where the SM radiative corrections are expressed in terms of the effective couplings 
Agl and As^ [0, and the top-quark mass dependence of the Zh^hj^ vertex 
correction in {g'i)sM is parametrized by the parameter 

- 175 GeV 

- 10 GeV • ^^-^^ 

The gauge boson propagator corrections, Agl and As^, are defined as the shift in 
the effective couplings gl{mlj and s^(m|J []T^ from their SM reference values at 
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= 175 GeV and rrijj = 100 GeV. They can be expressed in terms of the S and 
T parameters as 

Agl = ^|(m|J - 0.55635 = 0.00412Ar + 0.00005(1 - (100 GeV/m^:^)^], (3.6a) 
As^ = f{mlj - 0.23035 = 0.00360A5 - 0.00241AT - 0.00023a;„, (3.6b) 

where the expansion parameter is introduced to estimate the uncertainty of the 
hadronic contribution to the QED couphng l/a(mU = 128.75 ± 0.09 pT|] : 

l/a(m|J - 128.75 



0.09 



(3.7) 



Here, AS, AT, AU parameters are also measured from their SM reference values 
and they are given as the sum of the SM and the new physics contributions 



= A5sM + AT = ATsM + T,,,, AU = AUsu + U, 



The SM contributions can be parametrized as 



ASsM = -0.007Xi + 0.091x^ - 0.010x|^, 

ATgM = (0.130 - 0m3xH)xt + OmSxt - 0.079xj:^ - 0.028x^ 

+0.0026x1^, 
Af/gM = 0.022x4 - 0.002a;^, 



where x^j is defined by 



= log(m^/100 GeV). 



(3.8) 

(3.9a) 

(3.9b) 
(3.9c) 

(3.10) 



The pseudo-observables of the Z-pole experiments are given by using the above 
eight effective couplings as follows. The partial width of Zi boson is given by 



/ 



9l + 9r 



2 C 



fv 



+ 



f f 
9l~9r 



2 C 



fA 



1 + -Q, 



TT 



(3.11) 



where the factors Cfv and Cja account for the finite mass corrections and the final 
state QCD corrections for quarks. Their numerical values are listed in Table ^. 
The a^-dependence in Cqv, CqA is parametrized in terms of the parameter 

«s(mzj -0.118 



X, 



0.003 



(3.12) 



The last term proportional to a(m|J/7r in eq. ( p.ll| ) accounts for the final state 
QED correction. The total decay width Tz^ and the hadronic decay width Th are 
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u 


3.1166 + 0.0030a;, 


3.1351 + 0.0040X, 


d = s 


3.1166 + 0.0030a;, 


3.0981 + 0.0021X, 


c 


3.1167 + 0.0030a;, 


3.1343 + 0.0041X, 


b 


3.1185 + 0.0030a;, 


3.0776 + 0.0030X, 


V 


1 


1 


e = /i 


1 


1 


r 


1 


0.9977 



Table 4: Numerical values of factors Cfy^ ^fA for quarks and leptons used in eq. ( 3.11 ). 
The finite mass corrections and the final state QCD corrections for quarks are taken into 
accounted. 



given in terms of Fj: 

Tzi = 3r, + Te + + r, + r,,, 
= r„ + r, + + r, + r^. 

The ratios i?^, i?^, and the hadronic peak cross section cr° are given by: 



r,' ""^"mi n ■ 



(3.13a) 
(3.13b) 



(3.14) 



The left-right asymmetry parameter is also given in terms of the effective 
couplings gl^ as 

Af = (3.15) 
The forward-backward (FB) asymmetry A^p^^ and the left-right (LR) asymmetry 

f 

A£j^ are then given as follows: 



A^n = Af. 



(3.16a) 
(3.16b) 



3.2 W boson mass 

The theoretical prediction of can be parametrized as [|1^, |20[ 

m^(GeV) = 80.402 - 0.288 AS + 0.418 AT + 0.337 AU + 0.012 x^, (3.17) 



by using the same parameters, AS, AT, AU ( |3.8|) and (|3.7| ). 
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3.3 Observables in low-energy experiments 

In this subsection, we show the theoretical predictions for the electroweak observ- 
ables in the low-energy neutral current experiments (LENC) — (i) polarization 
asymmetry of the charged lepton scattering off nucleus target (|3.3.lH373l4D , (ii) 
parity violation in cesium atom ( p.3.5|) , (iii) inelastic z/^-scattering off nucleus tar- 
get ( p.3.6|) and (iv) neutrino-electron scattering ( p.3.7|) . The experimental data are 



summarized in Table Theoretical expressions for the observables of (i) and (ii) 
are conveniently given in terms of the model- independent parameters Cig, p2[ 
and ^ . The z/^-scattering data (iii) and (iv) are expressed in terms of the pa- 
rameters g'l'a ■ model-independent parameters can be expressed compactly 



in terms of the reduced hehcity amplitudes M^^ [|^, of the process faf'p faf'p'- 



Ci. = ^^( M[\ + M[\-M'^,-M'^^), (3.1a) 
C2, = ^^( M[\-Mi\ + M'^,-M'^^l (3.1b) 
Cs, = -^{-Mi\ + Ml^^ + - M^'^«), (3.1c) 



F 

gZ' = (3.id) 

Below, we divide these model-independent parameters into two pieces as 

Ciq = (Cjg)sM + '^Ciq, (3.2a) 
^?/ = teV + A^?/, (3.2b) 

where the first term in each equation is the SM contribution which is parametrized 
conveniently by AS* and AT in ref. [^. The terms AC^g and ^g'l'J represent the 
additional contributions from the Z-Z' mixing and the Z2 exchange: 



AC,^ = (-0.19., - 0.15c, + 0.655)^- 4 (^i_^I(%±^, (3.3a 



AC,a = (0.36s, - 0.54c, + 0.175)^- gl ^)(Og + Qg) ^ ^3 3^^ 



AC^^ = (0.02s, - 0.84c, + 1.485)^"- ^ „ ^^^/^ , , (3.3c) 



AC,, = (0.02s, + 0.84c, - 1.485)e - 7!^ ^ /^^ 2 ' ^^'^^^ 



cl l^plGpm^ 








cl 1^Gpm\ 




gl + 
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AC,^ = (-0.82c, + 1.005)e - f 

AC,, = (1.06., - 0.82c, - 1.00^)e - ^ 2V2G,ml ' 



^9 

^9[ 
^9 
^9 



LL — 



uu 



'LR 



ud 

LR — 



ue 

'LL = 



'lr 



0.44., + 0,22e,-0,m)C>|^^^ 
-0.35., + 0,01c, + 0,82^)f+U|L-. 
0.04.,- 0.72c, + 0.59^)f+|^-||ij-. 
-0.22.,- 0.52c,-0.4yK>U|ij-. 
0.12.,, + 0.28c, -0.23^)f+|^-|||j-. 
^(-0.14., + 0.49c,-1.23.)f+|^Jgi^. 



(3.3e) 
(3.3f) 
(3.3g) 
(3.3h) 
(3.3i) 
(3.3j) 
(3.3k) 
(3.31) 



where c, = cosPe and s, = 
3.3.1 SLAC eD experiment 

The parity asymmetry in the inelastic scattering of polarized electrons from the 
deuterium target was measured at SLAC The experiment constrains the 

parameters 2Ci„ — Cu and 2C2u — C2d- The most stringent constraint shown in 
Table H is found for the following combination 



AsLAc = 2Ci„ - + 0.206(2C2„ - C^d) 
= 0.745 - 0.016 Z\5 + 0.016 zAT 

+ 2ACi„ - ACi, + 0.206(2AC2„ - ACa^) 



(3.4a) 
(3.4b) 



where the theoretical prediction Q is evaluated at the mean momentum transfer 
(Q2) = 1.5 GeV2. 

3.3.2 CERN /i^C experiment 

The CERN fi^C experiment [^] measured the charge and polarization asymmetry 
of deep-inelastic muon scattering off the ^^C target. The mean momentum transfer 
of the experiment may be estimated at (Q^) = 50 GeV^ [p5|. The experiment 
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constrains the parameters 2C2u—C2d and 2C3„— Cs^. The most stringent constraint 
is found for the following combination |^ 

^CERN = 2C3„ - C3d + 0.777(2C2„ - C2d) (3.5a) 
= -1.42 - 0.016 AS + 0.0006 A T 

+ 2AC3„ - AC3rf + 0.777(2AC2„ - AC2,). (3.5b) 



3.3.3 Bates eC experiment 

The polarization asymmetry of the electron elastic scattering off the ^^C target 
was measured at Bates The experiment constrains the combination 



^Bates — Ciu + Cid (3.6a) 

= -0.1520 - 0.0023 AS + 0.0004 AT + ACi„ + ACid, (3.6b) 
where the theoretical prediction is evaluated at (Q^) = 0.0225 GeV^. 

3.3.4 Mainz eBe experiment 

The polarization asymmetry of electron quasi-elastic scattering off the ^Be target 
was measured at Mainz |^ . The data shown in Table |^ is for the combination 



^Mainz = -2.73Ci„ + 0.65Cid - 2. 19C2„ + 2.03C2d (3.7a) 
= -0.876 + 0.043A5 - 0.035Ar 

- 2.73ACi„ + 0.65ACirf - 2.19AC2„ + 2.03AC2d, (3.7b) 

where the theoretical prediction is evaluated at (Q^) = 0.2025 GeV^. 
3.3.5 Atomic Parity Violation 

The experimental results of parity violation in the atom are often given in terms of 
the weak charge Qy^{A, Z) of nuclei. By using the model-independent parameter 
Cj^g, the weak charge of a nuclei can be expressed as 

Q^(A, Z) = 2ZC,^ + 2{A - Z)C,^. (3.8) 



By taking account of the long-distance photonic correction ||2^, we find Cip and 

Cin as 

Cip = 0.03601 - 0.00681 A5 + 0.00477 AT + 2ACi„ + ACm, (3.9a) 
Ci„ = -0.49376 - 0.00366 AT + ACiu + 2ACm- (3.9b) 
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The data for cesium atom H Cs Q is given in Table ^ and the theoretical 



prediction of the weak charge is found to be 

QwilfCs) = -73.07 - 0.749 AS - 0.046 AT + 37QACiu + 422ACirf. (3.10) 

3.3.6 Neutrino-quark scattering 

For the z/^-quark scattering, the experimental results up to the year 1988 were 
summarized in ref. |^ in terms of the model-independent parameters gj^, 6l, 5\. 



The most stringent constraint on the result in ref. is found for the following 
combination: 

i^FH = 9l + 0.879^^ - 0.0105^ - 0.0434- (3.11) 

More recent CCFR experiment at Tevatron measured the following combina- 
tion [B2l 



KccFK = l-7897^i + 1.1479^ - 0.09165^ - 0.0782^. (3.12) 

The data are shown in Table ^ and the SM predictions are calculated from our 
reduced amplitudes ( |3.1cl| ) as follows 0, 



9i = i9zy+{9zr, si = {gzy-{9zr, (3.13) 

for a = L and R, respectively, where 

glf = 0.3468 - 0.0023AS + 0.0041AT, (3.14a) 

gl^^ = -0.1549 - 0.0023A5 + 0.0004AT, (3.14b) 

gl^^ = -0.4299 + 0.0012A5 - 0.0039Ar, (3.14c) 

gl^^ = 0.0775 + 0.0012A5 - 0.0002AT. (3.14d) 

The above predictions are obtained at the momentum transfer (Q^) = 35 GeV^ 
relevant for the CCFR experiment The estimations are found to be valid 



also for the data of ref. [^, whose typical scale is (Q^) = 20 GeV^. 
3.3.7 Neutrino-electron scattering 

The z/^-e scattering experiments measure the neutral currents in a purely leptonic 
channel. The combined results p, ^ are given in Table |^. The theoretical pre- 
dictions 

^^'f = -0.273 + 0.0033 AS - 0.0042 AT + A^^f , (3.15a) 
gl""^ = 0.233 + 0.0033 AS - 0.0006 AT + Ag^""^, (3.15b) 
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are evaluated at (Q^) = 2meE^ with = 25.7 GeV for the CHARM-II experi- 



ment 33 



4 Constraints on Z' bosons from electroweak ex- 
periments 

Following the parametrization presented in Sec. 3, we can immediately obtain 
the constraints on T^ew,C, and fi'|;/c^^l2 from the data listed in Table 2 and Ta- 
ble ^. Setting S'new = f^new = 0, wc find that the Z-pole measurements constrains 
Tnew and ^ while mw data constrains T^ew The contact term gE/c^^'zi 
strained from the LENC data. The number of the free parameters is, therefore, 
six: the above three parameters and the SM parameters, m^, a;s(mzi) and a(m|J. 
Throughout our analysis, we use 

= 175.6 ± 5.5 GeV [19], (4.1a) 
asiruzj = 0.118 ± 0.003 [|I3], (4.1b) 
l/a(m|J = 128.75 ± 0.09 [0, (4.1c) 



as constraints on the SM parameters. The Higgs mass dependence of the results 
are parametrized by x^j ( |3.10| ) in the range 77 GeV < <^ 150 GeV. The lower 



bound is obtained at the LEP experiment [Q. The upper bound is the theoretical 
limit on the lightest Higgs boson mass in any supersymmetric models that accom- 



modate perturbative unification of the gauge couplings We first obtain the 
constraints from the Z-pole experiments and W boson mass measurement only, 
and then obtain those by including the LENC experiments. 

4.1 Constraints from Z-pole and mw data 

Let us examine first the constraints from the Z-pole and mw data by performing 
the five-parameter fit for Tnew, ^, '"^t, «s("^Zi) and a(mU. The results for the 
Xi'4^iV and V models at 5 = are summarized as follows: 

(i) x-model (5 = 0) 

Tnew = -0.040 + 0.15x^± 0.12 1 ^ 

e = 0.00017 - 0.00005x^ ± 0.00046 J ^"^'^ ' (4.2) 

xLn/(d-o.f.) = (16.5 + 0.7x^)/(12), 
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(ii) ip-mode\ {6 = 0) 



T„ew = -0.043 + 0.16x;,± 0.11 1 ^ 

C = 0.00019 + 0.00012xj:^ ± 0.00050 J ' (4.3) 

xL„/(ci.o.f.) = (16.5 + 0.4x^)7(12), 

(iii) 77-model (S = 0) 

T,ew = -0.053 + 0.14x^± 0.11 1 ^ 

^ = -0.00014 - 0.00062x^ ± 0.00108 J ' (4.4) 

xLn/(d.o.f.) = (16.6 + 0.4x^)7(12), 

(iv) z/-model {6 = 0) 

T„ew = -0.042 + 0.15x^ + 0.11 1 ^ 

^ = 0.00016 + 0.00007x^ + 0.00042 j ' (4.5) 

xLn/(d.o.f.) = (16.5 + 0.5x^)7(12). 



In the above four Z' models, the resuhs for Tnew and ^ are consistent with zero for 
Xfj = 0. Moreover, the best fits of Tnew in all the Z' models are in the unphysical 
region, Tnew < 0. The parameter Tnew could be positive for the large Xj^: For 
example, Xjj = 0.41 {rrijj = 150 GeV) makes Tnew in all the four Z' models 
positive. The allowed range of the effective mixing angle ^ is order of 10"'^ for the 
?7-model and 10~^ for the other three models in 1-a level. The Xji^-dependence of 
^ in the ?7-model is larger than the other three models. For comparison, we show 
the result for the leptophobic //-model (5 = 173) 

(v) leptophobic r/- model {6 = 1/3) 

Tnew = -0.049 + 0.15x^ + 0.11 1 ^ 

^ = 0.00269 + 0.00026xj:^ + 0.00309 J ' (4.6) 

xL„7(d.o.f.) = (15.9 + 0.5x^)7(12). 

By comparing the 77- model with no kinetic mixing {S = 0) in eq. ( [4. 41) , we find 
significantly weaker constraint on ^. In Fig. |l], we show the 1-a and 90% CL 
allowed region on the (^, Tnew) plane in the r^-model with 6 = and 173 for 
mjj = 100 GeV. 

The best fit results at m^j = 100 GeV under the constraint Tnew > are shown 
in Table 2. We can see from Table 2 that there is no noticeable improvement of 
the fit for the V and u models at 5 = 0. The Xmm remains almost the same as 
that of the SM, even though each model has two new free parameters, Tnew and 
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Figure 1: The 1-a and 90% CL allowed region on the T^ew) plane in the r^-model 
with 6 = and 1/3. The shaded region (T^ew < 0) corresponds to the unphysical 
region. 

^. The fit slightly improves for the leptophobic r^-model (5 = 1/3) because of the 
smaller pull factor for the Rt data. The probability of the fit, 18.7% CL, is still 
less than that of the SM, 26.2% CL, because the Xmin reduces only 0.8 despite two 
additional free parameters. 

We explore the whole range of the parameters, Pe and 6. In Fig. H, we show 
the improvement in Xmin o^^r the SM value, Xmin(SM) = 16.9 (see Table 2): 

Ax'^xLnWE,S)-xLniSM), (4.7) 

where XminiPE, ^) is evaluated at the specific value of Pe and 6 for rrifj = 100 GeV. 
As we seen from Fig. 0, the Xmin depends very mildly in the whole range of the 
Pe and 6 plane, except near the leptophobic r^-model {Pe = tan~^(-y/5/3) and 
6 = 1/3) ||10[- Even for the best choice of Pe and 6, the improvement in Xmm 
is only 1.5 over the SM. Because each model has two additional parameters T^ew 
and ^, we can conclude that no Z' model in this framework improves the fit over 
the SM. The "x" marks plotted in Fig. ^ show the specific models which we will 
discuss in the next section. 
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Figure 2: Contour plot of Ax^ = xLniP^^) - Xmin(SM) for mjj = 100 GeV. The 
mixing angle /Se for the v models are shown by vertical dotted lines. 

The step of each contour is 0.2. The "x" marks on the plot show the specific 
models listed in Table || in Sec. 5. 

4.2 Constraints from Z-pole + niw + LENC data 

Next we find constraints on the contact term glE/c^Tn^^ by including the low- 
energy data in addition to the Z-pole and niw data. Because T^ew and ^ are 
already constrained severely by the Z-pole and mw data, only the contact terms 
proportional to glt/c^Tn^z^ contribute to the low-energy observables, except for the 
special case of the leptophobic //-model (5 = 1/3). 

We summarize the results of the six-parameter fit for the iIj,x,V and u models: 



(i) x-model 



T 



-0.063 +0.14x^ ±0.11 



^ = -0.00005 -0.00006x^± 0.00044 Vpc 



y(d.o.f.) 



0.26 



O.Olx 



H 



±0.21 



A mm/ 



;i9.9±0.9x^)/(20), 



'1.00 0.25 0.09\ 

1.00 0.15 ,(4.8a) 
1.00 J 

(4.8b) 



(ii) ■?/'-model 
T 

new 

gl/clm% 



-0.065 ±0.15a;^ ±0.11 
-0.00014 ± 0.00012a; ± 0.00050 ] 
1.66 ±0.19x^ ±2.90 



'1.00 0.19 0.07^ 
1.00 0.03 
1.00 



,(4.9a) 
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xLn/(d.o.f.) = (21.1 + 0.8x^)/(20), (4.9b) 
(iii) 77-model 



Tnew = -0.074 +0.14a;^ ±0.11 

^ = -0.00038 - 0.00063s^ ± 0.00106j>p^ 
^l/c>k=-0.62 +0.08a;^ ±0.87 



/1. 00 0.06 -0.05\ 

1.00 -0.22 ,(4.10c 

V 1.00 y 



xLn/(d.o.f.) = (20.8 ± 0.5x^)/(20), (4.10b) 
fiv) //-model 



Tnew =-0.061 ±0.15Xj;^ ±0.11 ^ 

^ = 0.00010 ±0.00006x^± 0.00041 
c,|/c2m|^ = -0.65 ±0.04xj;^ ±0.54 

xL„/(d.o.f.) = (20.1±0.8x^)/(20). 



1.00 0.21 0.07 
}p,,,,= \ 1.00 0.03 
1.00 




The contact term S^l^/c^^^l^ in the and rj models is consistent with zero in 1-cr 
level. Both the best fit and the l-cx error of the parameters T^ew and ^ in all the Z' 
models are slightly affected by including the LENC data: The best fit value of Tnew 
in all the Z' models cannot be positive even for the m^j = 150 GeV {x^j = 0.41). 
Since the leptophobic ?7-model does not have the contact term, the low-energy 
data constrain the same parameters Tnew and ^. After taking into account both 
the high-energy and low-energy data, we find 

(v) leptophobic r] model {S = 1/3) 

Tnew = -0.074 ±0.148x^± 0.110 1 

^ = 0.00157 ±0.00019xj:^± 0.00279 J ' (4.12) 

xLn/(d.o.f.) = (21.2 ± 1.0x^)/(21). 



The allowed range of ^ is slightly severe as compared to eq. 

The best fit results for rrifj = 100 GeV under the condition Tnew > are shown 
in Table ^. It is noticed that the best fit values for the weak charge of cesium atom 
ll^Cs in the x^V and z/ models are quite close to the experimental data. These 
models lead to Ax^ = —1.8 (x), —0.8 (r/) and —1.6 (z/). No other noticeable point 
is found in the table. 

The above constraints on g%/c^^fn\^ from the LENC data give the lower mass 
bound of the heavier mass eigenstate Z2 in the Z' models except for the leptophobic 
?7-model. In Fig. |^, the contour plot of the 95% CL lower mass limit of Z2 boson 
from the LENC experiments are shown on the (/?£,5) plane by setting qe = Qy 
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Figure 3: Contour plot of the 95% CL lower mass limit of the Z2 boson obtained 
from the LENC experiments for = Qy and 771^ = 100 GeV. The vertical dotted 
lines correspond to the ^ind v models. The limits are given in the unit of 

GeV. The "x" marks on the plot show the specific models listed in Table § in 
Sec. 5. 



and = 100 GeV under the condition m^a > 0. In practice, we obtain the 95% 
CL lower limit of the Z2 boson mass mgs in the following way: 

C95 dmz^P{mz^) 



0.05 



(4.13) 



where we assume that the probability density function Pijnz.^) is proportional to 
exp(-x2(m^J/2). 

We can read off from Fig. |^ that the lower mass bound of the Z2 boson in the 
if) model at 5 = is much weaker than those of the other Z' models. It has been 
pointed out that the most stringent constraint on the contact term is the APV 
measurement of cesium atom 0. Since all the SM matter fields in the ip model 
have the same U(l)' charge (see Table 1), the couplings of contact interactions are 
Parity conserving, which makes constraint from the APV measurement useless. 
We also find in Fig. |^ that the lower mass bound of the Z2 boson disappears near 
the leptophobic r^- model {Pe = tan" 



'5/3) and 6 = 1/3) [|Tg]. 
We summarize the 95% CL lower bound on mz2 for the ^^^d v models 

{5 = 0) in Table |^. For comparison, we also show the lower bound of in the 
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Table 5: The 95% CL lower bound of m^a (GeV) in the Xi'4'il ^-iid v models 
(5 = 0) for qe = Qy and rrifj = 100 GeV. The results of previous study |^ and 
of recent direct search are shown for comparison. 





X ^ 77 1/ 


Our results 
Langacker et al. [36 
direct search p7| 


451 136 317 284 
330 170 220 — 

595 590 620 — 



previous study [|36[ in the same table. The bounds on the and Zi, masses are 
more severely constrained as compared to ref. [36|. Although we used the latest 
electroweak data, our result for the Z^ boson mass is somewhat weaker than that 
of ref. [3^. In the analysis of ref. |3^, the e'^e" fi^fi^,T^T^ data below the Z 
pole are also used besides the Z-pole, mw and the LENC data. As we mentioned 
before, the lower mass bound of the Z' boson is obtained from the LENC data, 
not from the Z-pole data. Because the APV measurement which is most stringent 
constraint in the LENC data does not well constrain the ip model, it is expected 
that the e"'"e~ annihilation data below the Z-pole play an important role to obtain 
the bound of Z^ boson mass. 

Our results in Table ^ are also slightly weaker than those in ref. ^ . The results 
in ref. have been obtained without including the Z-Z' mixing effects and by 
setting = 175 GeV, = 100 GeV, = and Tnew = 0. 



4.3 Lower mass bound of Z2 boson 

We have found that the Z-pole, mw and the LENC data constrain (Tnew, 5 ^new 
and 5'|;/c^'^|2, respectively. We can see from eq. ( p. 17] ) that, for a given (, con- 
straints on Tncw5^ and glE/c^m^z^ can be interpreted as the bound on m^^. We 
show the 95% CL lower mass bound of the Z2 boson for 771^ = 100 GeV in four Z' 
models as a function of (. The bound is again found under the condition > 0. 
Results are shown in Fig. §.(a) ~ i[(d) for the V'? ^ models, respectively. The 
lower bound from the Z-pole and mw data, and that from the LENC data are 
separately plotted in the same figure. In order to see the (^s-dependence of the m^a 
bound explicitly, we show the lower mass bound for the combination mz2gY / 9e- 
We can read off from Fig. |^ that the bound on mz2gY/gE is approximately in- 
dependent of qe for Qe/qy = 0.5 ~ 2.0 in each model. As we expected from 
the formulae for Tnew and ^ in the small mixing limit (eq. ( |2.16| )), the Z2 mass 
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(a) X model 



(b) v model 




Figure 4: The 95% CL lower mass limit of Z2 in the x? ^) V v models for 
= 100 GeV. The Z2 boson mass is normalized by Qe/ Qy- Constraints from 
Z-pole experiments and LENC experiments are separately shown. The results of 
the direct search at Tevatron |37] for the Xi ^ ^"^^ V niodels are also shown. 

is unbounded from the Z-pole data at C = 0. For models with very small (, the 
lower bound on 'mz2, therefore, comes from the LENC experiments and the di- 
rect search experiment at Tevatron. For comparison, we plot the 95% CL lower 
bound on obtained from the direct search experiment |^ in Fig. ^. In the 
direct search experiment, the Z' decays into the exotic particles, e.g., the decays 
into the light right-handed neutrinos which are expected for some models, are not 
taken into account. We summarize the 95% CL lower bound on for the "4^, V 
and u models {6 = 0) obtained from the low-energy data and the direct search 
experiment in Table ||. 

The lower bound of is affected by the Higgs boson mass through the T 
parameter. As we seen from eqs. (|4.2|) ~ (^TSf), T^ew tends to be in the physical 
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region (Tnew > 0) for large m^j {x^). Then, we find that the large Higgs boson 
mass decreases the lower bound of m^^. For C = 1) the lower bound in the 
Xy'^^y^ iv) models for rrifj = 150 GeV is weaker than that for rrifj = 100 GeV 
about 7% (11%). On the other hand, the Higgs boson with m^j = 80 GeV makes 
the lower bound in all the Z' models severe about 5% as compared to the case 
for = 100 GeV. Because Tnew and ^ are proportional to ("^ and (, respectively 
(see eq. ( p.l6|) ), and it is unbounded at \(\ ~ 0, the lower bound of niay be 
independent of rrijj in the small \(\ region. The m^i^-dependence of the lower mass 
bound obtained from the LENC data is safely negligible. 

It should be noted that, at C = 0, only the leptophobic //-model (5 = 1/3) is not 
constrained from both the Z-pole and the low-energy data. The precise analysis 
and discussion for the lower mass bound of the Z2 boson in the leptophobic re- 



model can be found in ref. |^8[. It is shown in ref. |38| that the ^-dependence of 
the lower mass bound is slightly milder than that of the ?7-model with 5 = in 
Fig. |.(c). 

It has been discussed that the presence of Z2 boson whose mass is much heavier 
than the SM Z boson mass, say 1 TeV, may lead to a find-tuning problem to 
stabilize the electroweak scale against the U(l)' scale The Z2 boson with 

^ 1 TeV for qe = Qy allowed by the electroweak data only if ( satisfies the 
following condition: 

-0.6^C^+0.3 for the models, (4.} a) 

-0.7 ^ + 0.6 for the t] model. ^ ' ^ 

In principle, the parameter ( is calculable, together with the gauge coupling qe, 
once the particle spectrum of the Eq model is specified. In the next section, we 
calculate the ( parameter in several Eq Z' models. 



5 Light Z' boson in minimal SUSY ^^g-models 

It is known that the gauge couplings are not unified in the Eq models with three 
generations of 27. In order to guarantee the gauge coupling unification, a pair 
of weak-doublets, H' and H', should be added into the particle spectrum at the 
electroweak scale |^0[. They could be taken from 27 -|- 27 or the adjoint represen- 
tation 78. The U(l)' charges of the additional weak doublets should have the same 
magnitude and opposite sign, a and —a, to cancel the U(l)' anomaly. In addition, 
a pair of the complete SU(5) multiplet such as 5 + 5 can be added without spoiling 
the unification of the gauge couplings [0, 
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Table 6: Charge assignment for the extra weak doublets in the minimal Eq model 
and the ?7bkm niodel of ref. [^. The symbol ai{—ai) for i = x,ip,Vy^ the U(l)' 
charge of L or Hd {Hy). 





field 




1 






minimal model 


H' 


1 

2 




dip 




ay 




h' 










-ay 


r/BKM model [0 


H[ 


1 
2 






1 






Hi 


+ \ 






-1 








1 
2 






1 






h; 


+ i 






-1 






D' 


1 






2 






3 






3 








+1 






2 
3 





The minimal Eq model which have three generations of 27 and a pair 2 + 2 
depends in principle on the three cases; H' has the same quantum number as L or 
Hd of 27, or Hu of 27. All three cases will be studied below. The hypercharge 
and U(l)' charge of the extra weak doublets for the x, i^, ^ models are hsted in 
Table |^. For comparison, we also show those in the model of Babu et al. [p!0| , 
where two pairs of 2 + 2 from 78 and a pair of 3 + 3 from 27 + 27 are introduced 
to achieve the quasi-leptophobity at the weak scale. 

Let us recall the definition of 



c 



9z m^zz' 



(5.1) 



9e mz 

In the minimal model, the following eight scalar-doublets can develop VEV to give 
the mass terms m| and ml^' m eq. ( p.4| ): three generations of Hu,Hd, and an 
extra pair, H' and H'. Then, m| and m'^zz' written in terms of their VEVs as 



rrir 



:9l 



r 3 



m 



zz' —9z9e 



1=1 



(5.2a) 



(5.2b) 



where i is the generation index. The third component of the weak isospin 1^ for 
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the Higgs fields are 

hiH,) = h{H') = -h{H^) = -h(W) = 1/2. (5.3) 

Taking account of the U(l)' charges of the extra Higgs doublets, = —Qe', we 
find from eq. ( ^.l]) 

3 



C = ^- ^ - S. (5.4) 



i=i ^ ^ 



We note here that the observed /x-decay constant leads to the following sum rule 

1 



vl + vl + v% +vj^^v' = ^ (246 GeV)^ (5.5) 



where 



7)2 3 „,2 

i=i ^ i=i ^ 



(5.6) 



By further introducing the notation 



,2 I „,2 
~r 

V 

we can express eq. (|5.4|) as 



tan/3 = — , (5.7a) 

Vd 

= (5.7b) 



( = 2|-gf"(l -x2)sin2/? + g^'*(l-x2)cos2/? + gf'a;2| -5. (5.8) 

Because H' and if' are taken from 27 + 27, the U(l)' charge of H', Qe\ is 
identified with that of L, or Hu- 

Among all the models, only in the x-iiiodel one can have smaller number of 
matter particles. In the x-iiiodel, three generations of the matter fields 16 and a 
pair of Higgs doublets make the model anomaly free. In this case, ( is found to be 
independent of tan/?: 



C = 2 , , ,_2o ^ (5-9a, 



gg" -gfnan ^/j 
1 + tan^ (3 

2Ql' - 6. (5.9b) 
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Table 7: Coefficients of the 1-loop /3-functions for the gauge couphngs in the 
MSSM, the minimal Eq models and the ?7bkm niodel [|TU|. The model x(16) has 
three generations of 16 and a pair 2 + 2. The model x(27) and ip, Vy ^ have three 
generations of 27 and a pair 2 + 2. 





MSSM 


X(16) 


X(27) 




V 


V 


r/BKM H 


bi 


33 
5 


33 
5 


48 
5 


48 
5 


48 
5 


48 
5 


53 
5 




1 


1 


4 


4 


4 


4 


5 


h 


-3 


-3 














1 


bE 




6 + ^ 

^ 10 


9 + — 

^ 10 


9 + * 


9 + fa' 


9 + fa2 


77 

5 


biE 










-|a 

5 


-%a 

5 


16 

5 



Let us now examine the kinetic mixing parameter 5 in each model. The bound- 
ary condition of 5 at the GUT scale is 5 = 0. The non-zero kinetic mixing term 
can arise at low-energy scale through the following RGEs: 



d 


1 


—rC^i = 
dt 


2^ 


d 


1 


dt''' = 






1 






dt 


2^ 



feittj, (5.10a) 
[bE + 2biE5 + b^5')al (5.10b) 



'IE 



biS)ai, (5.10c) 



where i = 1, 2, 3 and t = In/x. We define ai and as 

_5g^ _hgl 
"^ = 3 4^' "^ = 3 4^- 

The coefficients of the /3-functions for cti, and 5 are: 

6i = ^Tr(F2), = ^Tr(Ql), h^E = ^Tr(FQij). (5.12) 

From eq. ( p.lOc|) , we can clearly see that the non-zero 5 is generated at the weak 
scale if biE 7^ holds. In Table we list 61, bE and biE in the minimal ^ 
models and the ?7bkm niodel [|10|- As explained above, the x(16) model has three 
generations of 16, and the x(27) model has three generations of 27. We can see 
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Table 8: Predictions for gE and 5 at /i = mz^ in the minimal models and the t/bkm 





i,e coupnng 


gy is fixed as gy 


— U.Ou. 


model 


a 


9e 


9e/9y 


6 


, /'I «N 


o 

o 


u.ooo 


n Q8Q 

u.yoy 


0.066 




-2 


0.361 


1.010 


-0.044 


-v-l'97'\ 


Q 
O 


u.ooo 


n Q8Q 


u.uuu 




-2 


361 


1 010 


-0.044 


V 


1 

X 


0.364 


1.020 






2 


0.356 


0.999 


056 




-2 


0.356 


0.999 


-0.056 


V 


1/6 


0.366 


1.025 


0.018 




-2/3 


0.351 


0.982 


-0.071 


V 




0.361 


1.010 


0.044 




-^3/8 


0.353 


0.989 


-0.066 


'7BKM 113 




0.308 


0.862 


0.286 



from Table |^ that the magnitudes of the differences hi — 62 and 62 — &3 are common 
among all the models including the minimal supersymmetric SM. This guarantees 
the gauge coupling unification at /i = mcuT — 10^^ GeV. 

It is straightforward to obtain gEif^Zi) and for each model. The ana- 

lytical solutions of eqs. ( ^.10a|) ^ ( ^.10c|) are as follows: 



1 



1 1 , , rncuT 
1 bi in . 

acuT 27r mz^ 



>IE 



hi 



1 



+ 



ai[mz^, 



O-GUT 



bE 
hi 



biE 
hi 



1 



acuT L Oi \ Oi J ) laiiniz,] 
'hiEy aiimzi) - acuT 
hi J a ' 



acvT 



(5.13a) 
(5.13b) 

(5.13c) 



GUT 



where acuT denotes the unified gauge coupling at /x = mcuT- In our calculation, 
Oizijnzi) = 0.118 and a(m^J = e^(m^J/47r = 1/128 are used as example. These 
numbers give gy{mz^) = 0.357. We summarize the predictions for gE and 6 at 
fi = rrizi in the minimal Eq models and the ?7bkm model in Table |^. In all the 
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Table 9: Predictions for the effective Z-Z' mixing parameter in the minimal 
X, V and the ?7bkm niodel for = and 0.5, and tan /? = 2 and 30. 









= 


= 


= 0.5 






tan/5 


tan (3 




a 


2 


30 


2 


30 


X 


3 


-0 


.88 


0.14 




-2 




-0.77 






1 


0.60 


1.02 


0.55 


0.76 




2 


0.58 


1.00 


0.79 


1.00 




-2 


0.69 


1.11 


-0.16 


0.06 


V 


1/6 


-1.02 


-1.35 


-0.35 


-0.52 




-2/3 


-0.93 


-1.26 


-1.11 


-1.26 


V 




0.36 


0.77 


0.57 


0.77 




-a/3/8 


0.47 


0.88 


-0.34 


-0.14 






-1.29 


-1.62 


-1.79 


-1.95 



minimal models, the ratio QE/gv is approximately unity and \6\ is smaller than 
about 0.07. On the other hand, the ?7bkm model predicts gE/gr ~ 0.86 and 
6 ~ 0.29, which is close to the leptophobic-r] model at 5 = 1/3. In Figs. ^ and ^, 
we show the predictions of all the models by "x" symbol. 

Next, we estimate the parameter ( for several sets of tan/3 and x. In Table 
we show the predictions for ( in the minimal 4'^ Vy ^ models and the t^bkm model. 
The results are shown for tan (3 = 2 and 30, and = and 0.5. We find from the 
table that the parameter C, is in the range |C| ^ 1-35 for all the models except for 
the %KM model, where the predicted C lies between —2.0 and —1.2. It is shown in 
Fig. H that mz^gy/gE is approximately independent of gE/gv- Actually, we find 
in Table ^ and Table ^ that the predicted \S\ is smaller than about 0.1 and gs/gY 
is quite close to unity in all the minimal models. We can, therefore, read off from 



Fig. ^ the lower bound of m^a in the minimal models aX qe = gv- In Table |T0 
we summarize the 95% CL lower mz^ bound for the minimal Xi Vi ^ models and 
the r/BKM model which correspond to the predicted C in Table |[ Most of the 
lower mass bounds in Table [T0| exceed 1 TeV. The Z2 boson with uiz^ ~ 0(1 TeV) 
should be explored at the future collider such as LHC The discovery limit of the 
Z' boson in the models at LHC is expected as |^| 



3040 2910 2980 *** 
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Table 10: Summary of the 95% CL lower bound of (GeV) which corresponds 
to the predicted ( in Table ||. 
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30 
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30 
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1330 


620 




-2 




1230 
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1800 
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1480 
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1750 
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1750 




-2 


1380 


1890 


520 


370 




+ 1/6 


1330 


1690 


620 


790 




-2/3 


1230 


1590 


1410 


1590 


V 




1180 


1720 


800 


520 






1030 


1580 


1320 


1580 






1520 


1930 


2150 


2360 



All the lower bounds of niz^ listed in Table [1^ are smaller than 2 TeV and they 
are, therefore, in the detectable range of LHC. But, it should be noticed that most 
of them (1 TeV ^m^j) may require the fine-tuning to stabilize the electroweak 
scale against the U(l)' scale [0^] . 

The lower bound of the Z2 boson mass in the r^BKM model for the predicted ( 
can be read off from Fig. 2 in ref. |^8[. Because somewhat large ( is predicted in 
the r/BKM model, 1 ^ |C|, the lower mass bound is also large as compared to the 
minimal models. 



6 Summary 

We have studied constraints on Z' bosons in the SUSY Eq models. Four Z' models 
— the X; V u models are studied in detail. The presence of the Z' boson 
affects the electroweak processes through the effective Z-Z' mass mixing angle ^, 
a tree level contribution Tnew which is a positive definite quantity, and the con- 
tact term QElc^Tn^z^. The Z-pole, niw and LENC data constrain (rnew,^) ^new 
and 9% / (^\^'z2i respectively. The convenient parametrization of the electroweak 
observables in the SM and the Z' models are presented. From the updated elec- 
troweak data, we find that the Z' models never give the significant improvement 
of the x^"fit even if the kinetic mixing is taken into accounted. The 95% CL lower 
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mass bound of the heavier mass eigenstate Z2 is given as a function of the ef- 
fective Z-Z' mixing parameter C^. The approximate scahng low is found for the 
S's/S'y-dependence of the lower limit of rnz^- By assuming the minimal particle 
content of the model, we have found the theoretical predictions for C,. We have 
shown that the models with minimal particle content which is consistent with 
the gauge coupling unification predict the non-zero kinetic mixing term 5 and the 
effective mixing parameter of order one. The present electroweak experiments 
lead to the lower mass bound of order 1 TeV or larger for those models. 
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